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EFFECTS OF LIQUID DEPTH ON LATERAL SLOSHING 
UNDER WEIGHTLESS CONDITIONS 
by Jack A. Salzman, Thom A. Coney, and William J. Masica 
Lewis Research Center 

SUMMARY 

As a part of the general study of liquid behavior in weightlessness, an experimental 
drop-tower investigation was conducted to determine the effects of shallow liquid depths 
on lateral (asymmetric) sloshing in a zero Bond number environment. Contact angles 
were restricted to very near 0°, so that the sloshing equilibrium interface was hemi- 
spherical. Cylindrical tanks, 0. 635 and 0. 952 centimeter in radius, were tested with 
both hemispherical and flat bottoms. Natural frequency data were correlated as a func- 
tion of normalized liquid depth for both tank bottom shapes, and damping data were com- 
pared with predictions from high Bond number theory. 


INTRODUCTION 

Design criteria for space- vehicle stability and control include the ability to predict 
the behavior and characteristics of liquid-propellant systems when they are experiencing 
lateral sloshing motion. The frequencies and damping of the sloshing liquid mass must 
be understood under all gravity conditions and for a wide range of fillings (including very 
shallow liquid depths) in order to permit an analysis of possible interactions with ve- 
hicle stability and other systems of the vehicle. A comprehensive summary of most 
of the high Bond number slosh investigations and some of the original work by Reynolds 
and Satterlee in low Bond number sloshing is contained in reference 1. (Bond number is 
the ratio of acceleration to capillary forces. ) For high Bond numbers (greater than ap- 
proximately 50), where the equilibrium liquid surface is reasonably flat, analytical and 
experimental correlations exist for a wide range of sloshing phenomena pertinent to 
space-vehicle applications. For low Bond numbers, and especially with low contact 
angles where the equilibrium interface is highly curved, only limited information is 
available. 


Current information on low Bond number sloshing may be found in references 2 to 5. 
In general, these studies have considered certain basic characteristics of both free and 
forced lateral sloshing motion in low Bond number environments. With the exception of 
the analytical study in reference 3 and limited experimental work, no other investiga- 
tions of the effect of shallow liquid depths on low Bond number sloshing have been 
published. 

The purpose of this report is to present the results of an experimental investigation 
conducted at the NASA Lewis Research Center of the effects of shallow liquid depths on 
lateral sloshing under effectively zero (i. e. , less than 0. 001) Bond number conditions. 
Cylindrical tanks were used with both hemispherical and flat bottoms. Contact angles 
were restricted to very near 0°, so that the sloshing equilibrium liquid surface was 
nearly hemispherical. Further restrictions were that the slosh oscillations were of 
small amplitude and that the contact angle did not vary during sloshing. The data ob- 
tained on natural frequency are correlated as a function of normalized liquid depth by 
using the analytical results of reference 3. Damping data are compared with high Bond 
number correlations. 


SYMBOLS 


2 

a system acceleration, cm/sec 

B Q Bond number (B Q = aR 2 /j3) 

o 

g acceleration due to gravity, cm/sec 

h centerline interface height from tank bottom, cm 

K^, Kp Kg nondimensional constants 

R cylinder radius, cm 

t time, sec 

tj formation time, sec 

3 

V L liquid volume, cm 

3 2 

0 specific surface tension (j3 = a/p), cm /sec 

6 logarithmic decrement 

6 contact angle, deg 

p viscosity, cP 

O 

p density, g/ cm 
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a surface tension, dynes/cm 

o n o o 

SI natural frequency parameter or fundamental mode eigenvalue (fi = //3) 

natural frequency, rad/sec 


BACKGROUND INFORMATION 
Natural Frequency 

Flat-bottom cylinders. - If the liquid depth in a lateral (asymmetric) sloshing sys- 
tem is made shallow enough, the natural frequency of free oscillation will decrease. 

For the case of a flat-bottomed cylinder and a flat equilibrium liquid surface (i. e. , 90° 
contact-angle boundary condition), the free natural frequency may be expressed 
nondimensionally in the following form (ref. 1, p. 415): 

ft 2 = (6.2 55 + 1. 841 B q ) tanh ^1. 841 iij (1) 


where SI 


2 


is defined as 


SI 


2 


o> 2 R 3 


0 


( 2 ) 


and h is the centerline depth between the liquid surface and the tank bottom. For large 
Bond numbers, these equations reduce to the normal-gravity equation 

w 2 = 1. 841 £ tanh (l. 841 (3) 

R \ R/ 


No simple relation exists for other contact angles at arbitrary Bond numbers, par- 
ticularly for those combinations of contact angles and Bond numbers resulting in a highly 

p 

curved equilibrium liquid surface. At deep liquid depths (h/R > 2), the value of SI 
was experimentally determined to be 2. 6 for a zero Bond number, 0° contact- angle sys- 
tem (ref. 4). Limited experimental data for low- contact-angle, low Bond number sys- 
tems (ref. 6) indicate that the frequency dependence on liquid depth departs considerably 
from the form of equation (1). 

Hemispherical-bottom cylinders . - A recent analytical study of low Bond number 
(including zero) lateral sloshing in cylinders with hemispherical bottoms is contained in 
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Figure 1. - Computed, zero Bond number eigenvalues as function of 
liquid depth for cylindrical tank with hemispherical bottom. Contact 
angle, 5°; computed values from reference 3. 


reference 3. By using finite-difference techniques, normal modes were computed over 
a range of liquid depths at a contact angle of 5°. The computed, fundamental- mode, 
zero Bond number eigenvalues (ref. 3, p. 94) are plotted in figure 1 as a function of 
normalized liquid depth. As before, h is the centerline depth between the liquid surface 
and the tank bottom. A hyperbolic tangent approximation can be used to fit a functional 
relation to these calculated points. This relation, represented by the curve in figure 1, 
is 

O 2 = 2. 8 tanh ^2 (4) 

or, in terms of the natural frequency, 

to 2 = 2 . 8 JL tanh 1 2 — ^ (5) 

ft, 2 ' H / 

The empirical value of 2 in the argument of the hyperbolic tangent function is, of course, 
an approximation accurate to the indicated number of significant figures. It is interest- 
ing to note the similarity between the depth-dependent terms in equations (1) and (4). 

For depth ratios larger than about 2, bottom effects are negligible. The calculated 
zero Bond number eigenvalue of 2.8 (5° contact angle) at these deep depths agrees fa- 
vorably with the experimental value of 2. 6 (0° contact angle) reported in reference 4. 
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In this investigation the zero Bond number natural frequency dependence on liquid 
depth for both hemispherical-bottom and flat-bottom cylinders was obtained experimen- 
tally. The study was restricted to near 0° static contact angles that did not vary during 
the sloshing motion. Under these conditions, the equilibrium liquid-vapor interface was 
hemispherical in all cases. The results in hemispherical-bottom cylinders were ob- 
tained first and compared with the analytical study of reference 3 (i. e. , eq. (5)). An 
analogous equation describing the results for flat-bottom cylinders was then obtained 
empirically. 


Damping 

Deep liquid depth (h/R > 2) lateral slosh damping can be described in terms of the 
logarithmic decrement by 
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( 6 ) 


The value of is 6. 1 for high Bond numbers (ref. 7) and 28. 1 for zero Bond numbers 
(ref. 4). 

If the liquid depth is made sufficiently shallow, the damping is expected to increase. 
For high Bond numbers and flat-bottom cylinders, the increase in damping is predicted 
analytically in reference 8 and verified experimentally in reference 6. From these 
studies, the logarithmic decrement can be predicted by the following relation (ref. 1, 
p. 110): 


6 = 


5.23 M 1/2 p- 1/2 R- 3/4 g- 1/4 

tanh" 1 / 4 (l. 84 h.) 

1 + 2 / 

1 - — \csch / 

3. 68 — ^ 


\ R/_ 

V 

R/ \ 

R/_ 


( 7 ) 


Substituting the high Bond number natural frequency (eq. (3)) into equation (7) yields 



This equation shows that the total increase in high Bond number damping at shallow 
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depths is due to two factors: (1) the decrease in natural frequency to ^ contained in the 
first bracket of equation (8); and (2) the increased viscous effects caused by the proxim- 
ity of the tank bottom. The functional increase in damping due to the increased viscous 
effects is expressed by the second bracketed term in equation (8). 

High Bond number, shallow-depth damping relations for other tank geometries are 
presented in reference 1. In general, the same separation in factors as given in equa- 
tion (8) may be obtained, although the forms become more complex because of their em- 
pirical nature. Analytical or experimental investigations of the damping depth depend- 
ence in low Bond number environments appear nonexistent. 

If the logarithmic decrement is normalized in the form 



then the right side of the equation represents the increase in damping caused by viscous 
bottom effects alone. Thus, a direct comparison can be made between high Bond number 
and zero Bond number damping forms by appropriate substitution of the respective val- 
ues for K^. While this approach is not rigorously valid (in one case the interface is flat 
and in the other it is hemispherical), it is, nevertheless, useful in establishing compar- 
ative damping magnitudes. 


APPARATUS AND PROCEDURE 

The experimental apparatus used in this study is shown in figure 2. A near zero 
Bond number environment was obtained by permitting the experiment package to free 
fall a distance of 26 meters (a free fall time of approximately 2. 2 sec). Precautions 
taken to reduce the air drag on the package resulted in an effective acceleration level of 

5 

less than 10 g. This low level of acceleration together with the range of specific sur- 

o n 

face tensions (11. 8 to 28. 5 cm /sec ) and small cylinder radii (0. 635 and 0. 952 cm) 
used in the experiments resulted in Bond numbers of less than 0. 001. This was re- 
garded as a zero Bond number environment. 

Flat-bottom cylinders were assembled from precision-bore, borosilieate glass 
tubing and plate glass. Bottom edge distortions resulting from assembly were insignif- 
icant. Hemispherical- bottom cylinders were machined from acrylic plastic. The 
analytic- reagent-grade test liquids exhibited near 0° static contact angles on the cylinder 
surfaces. Care was taken to avoid a dynamic contact angle variation during free slosh- 
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TABLE I. - SUMMARY OF PARAMETERS 


Liquid Surface Density, a Viscosity, a Cylinder Normalized Measured Logarithmic 
tension, a p, p, radius, liquid natural decrement, 

a, g/ cm^ cP R, depth, frequency, 6 

dynes/ cm cm h/R u 

rad/sec 


Hemispherical-bottom cylinders 


Trichloro- 

18.6 

1. 58 

0. 70 

0. 635 

2.00 

10.6 

0. 99 

trifluoro- 





1.75 

11.0 

— 

ethane 





1. 50 

10.6 

— 






1.25 

10.3 

— 






1. 10 

10.8 

.75 






1. 00 

10.4 

1. 11 






.90 

10.8 

.87 






.75 

10.2 

.98 






.60 

9.23 

.84 






. 50 

9.34 

.96 






.40 

7.66 

1. 18 






.30 

7.48 

1.45 






.20 

6.96 

1. 47 

Carbon 

26. 9 

1. 59 

0.97 . 

0. 635 

2.00 

13.2 

0.85 

tetra- 





.85 

12.4 

1. 07 

chloride 





.65 

11.9 

.86 






.40 

10. 5 

.96 






.20 

8.26 

1.20 

Methanol 

22. 6 

0.79 

0.60 

0. 952 

2. 00 

9.40 

0.78 






.75 

8.84 

.77 






. 50 

7.95 

1.00 






.35 

7.22 

1. 06 






.25 

6.94 

1.47 


Flat- bottom cylinders 



L At 20° C. 


. 80 
.77 
1. 12 



ing, and for the data presented, no variation was observed. The tank radii and test liq- 
uids were chosen so as to maximize the number of sloshing oscillations while retaining 
a large enough cylinder size to ensure that the interface motion could be measured eas- 
ily. Liquid depths were calculated before each test by means of a volumetric relation, 
and these volumes were placed in the cylinders with a hypodermic syringe. Liquid 
depths were also measured directly during the data analysis. Properties of the test liq- 
uids, the cylinder sizes, and the liquid depths used are presented in table I. 

Small-amplitude, lateral sloshing motion was induced by initially displacing the 
liquid-vapor interface from its equilibrium position and then allowing its free motion. 
The lateral impulse which initially displaced the interface was supplied by an air piston 
(fig. 2) and was activated after a predetermined time increment which allowed the inter- 
face to form its zero Bond number configuration. The minimum time required for the 
interface to reach a sufficiently quiescent configuration was empirically determined to 
be (ref. 4) 



( 10 ) 


While this allowed formation time was not sufficient to ensure a completely quiescent 


Equilibrium interface 

Sloshing interface 




Equilibrium profile 


Sloshing profile 

(a) Hemispherical -bottom cylinder. 




Equilibrium profile 



Sloshing profile 
(b) Flat -bottom cylinder. 


Figure 3. - Liquid-vapor interface configurations. 
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Interface edge position relative to tank bottom, cm 


O 


2 . 200 ' 
2. 175 
2. 150 
2. 125 
2. 100 
2.075 
2.050 
2.025 
2.000 
1.975 



.8 .9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 

Time from initiation of free fall, sec 

Figure 4. - Sample data plot of interface edge motion. 



hemispherical interface, the motion was sufficiently damped so that it had no effect on 
the ensuing slosh wave. 

Interface motion was recorded by means of high-speed photography. Time measure- 
ments were obtained by viewing a precision digital clock with a calibrated accuracy of 
0. 01 second. Figure 3 shows both the equilibrium configuration and the lateral- sloshing 
(displaced) profile for both tank geometries tested. The initial process in obtaining the 
natural frequency and damping data involved the use of a motion- picture analyzer to re- 
duce the sloshing motion from pictorial to numerical form. Motion of the interface lead- 
ing edge at the tank wall was recorded as a function of time. The result, from which the 
data were obtained, is a graph (fig. 4) showing a damped sinusoidal wave for each edge. 


PRESENTATION OF DATA AND DISCUSSION OF RESULTS 
Natural Frequency 

Hemispheric al-bottom cylinders . - The first part of this investigation was directed 
toward experimental measurements of the zero Bond number natural frequency at shal- 
low liquid depths in hemispherical-bottom cylinders. Two different sized cylinders and 
three different liquids were tested, with liquid depth as the primary experimental pa- 
rameter. Data were obtained at liquid depths down to a minimum of 0. 2 radius. At 
depths below this minimum, the liquid-vapor interface was merely a thin film along the 
tank bottom, and measurements of the sloshing motion were not possible. The particu- 
lar systems which were tested and their respective frequency data are presented in 
table I. The natural frequency for each system (i. e. , for each set of values for R, j3, 
and h/R) was calculated by averaging the half periods for that test run. The recorded 
frequencies for this bottom shape had an average mean deviation of 6 percent. 

The natural frequency data are presented in figure 5 in normalized form so as to 
utilize the empirically determined, deep-liquid depth constant of 2. 6 (ref. 4). The curve 
in figure 5 is the approximated analytical depth dependence, tanh [2 (h/R)] , from equa- 
tion (5); the data agree favorably with this curve. Therefore, the zero Bond number 
natural frequency in hemispherical-bottom cylinders can be described by the relation 



A significant decrease in the natural frequency does not occur until the liquid depth is 
lowered below 1 radius. Note the similarity between the frequency-depth dependence of 
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f 



Figure 5. - Effect of liquid depth on natural frequency in hemispherical -bottom cylinders. 


equation (11) and of the flat- interface, flat-bottom case previously predicted in equa- 
tion (3). 

Although the frequency measured was not truly the natural frequency but was, in- 
stead, a damped-system frequency, the measured damping was sufficiently low so that 
the variation from the natural frequency was in all cases less than 5 percent. 

Flat-bot tom cylinders. - For the investigation of the effects of shallow liquid depths 
on the zero Bond number natural frequency in flat-bottomed cylinders, the same tank 
radii and test liquids were used as those for the hemispherical- bottom investigation. 
With this flat-bottom shape, data for liquid depths down to essentially zero were ob- 
tained. Measured frequencies for this bottom shape are recorded in table I. The re- 
corded frequencies had an average mean deviation of 5 percent. 

The frequency data are presented in normalized form in figure 6. The solid curve 
in this figure represents equation (11). The data in this figure indicate that the natural 
frequency for a flat- bottom cylinder starts to decrease at a lower liquid depth than does 
the natural frequency for a hemispherical-bottom cylinder and that the correlation for 
hemispherical-bottom cylinders (eq. (11)) cannot be used to represent the natural fre- 
quency for a flat-bottom cylinder. The functional trend shown by these data does indi- 
cate that a hyperbolic tangent relation can be used to represent the data if an additional 
constant is used. The resulting relation is of the form 


« 1 cc tanh 1 / 2 K x 


( 12 ) 
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If the value of Kj is taken to be 2 in order to facilitate comparison with the previous 
correlation, the value of K 2 best representing the data is 0.25. The dashed curve 
shown in figure 6 is the relation 



(13) 


Examination of the figure shows that the zero Bond number natural frequency in flat- 
bottom cylinders can be satisfactorily predicted by equation (13). Equation (13) and fig- 
ure 6 show that the natural frequency in flat-bottom cylinders does not decrease signifi- 
cantly until the liquid depth is lowered below 0. 75 radius. 

Equivalent bottom shapes. - High Bond number sloshing investigations (refs. 1 
and 9) have shown that tanks with different bottom shapes may be represented by a single 
equivalent bottom shape through normalization or equating of liquid volumes. The form 
of equation (12) and the relative magnitude of Kg (i. e. , 0. 25) suggest that the same ap- 
proach can be used for a normalization of the zero Bond number frequency data. Data 
obtained in the flat- bottom cylinders can be represented in equivalent hemispherical- 
bottom form by applying a value of 0. 33 to Kg. This equalizes the volumes for both 
bottom shapes at any value of h/R. Equation (12) is then in the form 


cc tanh 


1/2 


k 4 +# - 33 ) 


(14) 
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Figure 7. - Comparison of bottom shape effect on natural frequency by use of equivalent volumes. 


or in terms of normalized liquid volume 


w. a tanh 


1/2 



(15) 


Normalized frequency data from both bottom shapes are presented as a function of nor- 
malized liquid volume in figure 7. Although some parametric spread is evident, the 
correlation is sufficient for engineering purposes if a value of 2 is employed for K^. 

Equations (14) and (15) should not, however, be used to predict a sloshing frequency in 

0 

flat-bottom cylinders at a liquid volume V, less than 0. 33/7 tR (or h/R < 0) because 
this volume condition would obviously require a distortion of the hemispherical interface 
and of the sloshing mode. 


Damping Data 

Data on the effects of shallow liquid depths on zero Bond number slosh damping were 
obtained in both hemispherical- and flat-bottom cylinders. The damping was exponential 
for all tests and was measured in terms of the logarithmic decrement. Values were de- 
termined graphically from measurements of the decaying amplitudes. An example of 
this method may be found in reference 4. The logarithmic decrement for each set of 
parameters is given in table I. Because of the limited number of measured amplitudes 
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Zero Bond number data 



Figure 8. - Effect of liquid depth on damping. 


in each test run, there was considerable scatter in the damping data which prohibited 
any specific correlation. 

The damping data are shown in figure 8. These data have been normalized as in 
equation (9) by using the zero Bond number value of 28. 1 for and the actual meas- 
ured frequencies given in table I. As has been discussed in the section titled Damping, 
this was done so that the observed increase in damping can be attributed solely to vis- 
cous bottom effects and not to a decrease in the natural frequency. The curve plotted in 
figure 8 represents the correlated, high Bond number (flat bottom), viscous-damping 
depth dependence. This depth dependence has been normalized as in equation (9) by 
using 6. 1 for the value of and represents the second bracketed term in equation (8). 

Even though the data contain considerable scatter, it is evident from figure 8 that 
the depth dependence of zero Bond number damping cannot be correlated by the verified 
high Bond number theory. Although the damping does increase with decreasing liquid 
depth, a significant increase does not occur until a lower depth than that indicated by the 
high Bond number theory (i. e. , at an h/R of approximately 0. 5 rather than 1. 0). 

No damping data were obtained at a liquid depth of zero because of a very high re- 
duction in the slosh amplitude in only one cycle of oscillation at this depth. 


SUMMARY OF RESULTS 

An experimental investigation was conducted to determine the effect of shallow liq- 
uid depths (i. e. , less than 2 tank radii) on the natural frequency and damping character- 
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istics of lateral (asymmetric) sloshing in a zero Bond number environment. The study 
employed right circular cylinders, 0. 635 and 0. 952 centimeter in radius R, with both 
hemispherical and flat bottoms. Test liquids were restricted to those which possess 
near 0° static contact angles on the cylinder walls; the contact angle did not vary during 
the sloshing motion. Specific surface tensions |3 of the test liquids ranged from 11. 8 
to 28. 5 cubic centimeters per second squared. The study was restricted to small am- 
plitude oscillations and yielded the following results; 

1. The natural frequency in hemispherical-bottom cylinders was correlated as 
a function of normalized liquid depth. Experimental results can be described by the re- 
lation 



where h is the centerline height from the liquid-vapor interface to the tank bottom. 

2. The natural frequency in flat- bottom cylinders follows a relation similar to that 
in hemispherical-bottom cylinders. Experimental results for this bottom shape can be 
described by the relation 



3. The damping in terms of the logarithmic decrement increases with decreasing 
liquid depth for both bottom shapes. The depth at which damping starts to increase in a 
zero Bond number environment is lower than that at which a similar increase occurs at 
large Bond numbers. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 22, 1967, 
124-09-03-01-22. 
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